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PLASMA  INSTABILITIES  IN  THE  HIGH  AND  LOW  LATITUDE 
E-REGION  INDUCED  BY  HIGH  POWER  RADIO  WAVES 


1 .  INTRODUCTION 

The  suggestion  of  using  a  ground-launched  powerful  high-frequency  (HF) 
radio  wave  to  influence  (suppress  or  excite)  naturally  occurring  plasma 
instabilities  in  the  ionosphere  was  originally  made  by  Lee  et  al.  [1972] 
for  the  equatorial  electrojet  case.  Similar  work  had  been  carried  out 
later  by  Bujarbarua  and  Sen  [1978]  for  the  equatorial  F-region,  Stenflo 
[1983]  for  the  high-latitude  E-region  modes,  Keskinen  et  al.  [1983]  for 
convective  modes  in  the  high-latitude  F-region,  and  Chaturvedi  et  al. 
[1987a]  for  the  electrostatic  ion-cyclotron  (EIC)  instability  in  the  high- 
latitude  ionosphere.  These  papers  considered  the  parametric  excitation  of 
naturally  occurring  plasma  modes  by  the  high-frequency  pump  wave  with  the 
mode-coupling  arising  due  to  ponderomotive  force  (pmf)  effects.  For  the 
collisional  plasmas  (such  as  the  weakly  ionized  plasma  of  the  ionosphere), 
it  has  been  pointed  out  by  many  authors  that  the  nonlinear  coupling  via 
Joule  heating  of  electrons  is  an  important  effect  and  is  likely  to 
dominate  over  the  ponderomotive  force  effects  [Perkins,  1974;  Fejer,  1979; 
Gurevich,  1978;  Papadopoulos  et  al.,  1983].  This  latter  approach  has 
recently  been  incorporated  in  many  studies  such  as  those  by  Stenflo  [1985] 
for  convective  modes  at  high-latitudes,  Laxmi  and  Tripathi  [1987]  for 
equatorial  interchange  modes,  and  Sharma  [1989]  for  the  equatorial 
gradient-drift  modes,  among  others.  Chaturvedi  and  Ossakow  [in 
preparation,  1989]  have  also  included  these  effects  in  the  analysis  of 
convective  fluid  modes  in  the  high-latitude  ionosphere.  In  this  report, 
we  consider  the  effects  of  an  HF  pump  on  the  E-region  plasma  modes, 
applicable  to  the  high-  and  low-latitudes. 

There  have  been  many  experiments  using  an  HF  pump  to  study  its  effect 
on  the  high-latitude  E-region  [Stubbe  and  Kopka,  1983;  Hibberd  et  al., 
1983;  Coster  et  al.,  1985;  Djuth  et  al.,  1985;  Hoeg  et  al.,  1986;  Schlegel 
et  al.,  1987;  Wagner  et  al.,  1990,  among  others).  The  observations 
indicate  that  a  number  of  physical  effects  probably  occur  at  different 
times  (depending  apparently  on  ambient  conditions  and  also  on  the  pump 
wave  characteristics).  On  several  occasions,  in  the  above  mentioned 
experiments,  the  observed  effects  (primarily  irregularity  observations) 
could  be  attributed  to  a  presence  and/or  generation  of  naturally  occurring 
plasma  instabilities  in  the  high-latitude  E-region. 

Manuscript  approved  July  23,  1990. 


1 


In  this  paper,  we  extend  the  equatorial  electrojet  work  of  Lee  et  al. 
[1972]  to  the  high-latitude  E-region  case.  As  mentioned  earlier,  the 
nonlinear  Joule  heating  of  electrons  is  the  dominant  coupling  mechanism  in 
most  situations.  An  important  difference  in  the  study  of  this  region  from 
the  equatorial  case  is  the  possibility  of  a  variety  of  naturally  occurring 
plasma  modes  possible  at  the  auroral  latitudes.  In  addition  to  the  two- 
stream  and  gradient-drift  instabilities  associated  with  the  electrojet 
(lower  E-region)  and  common  to  the  low  and  high-latitudes,  at  high- 
latitudes,  modes  driven  by  parallel  currents  (at  upper  E-region  altitudes) 
(ion  acoustic,  ion-cyclotron  and  current  convective  instabilities)  should 
also  be  considered.  Thus,  the  observed  effects  in  an  artificial 
ionospheric  modification  (AIM)  experiment  by  HF  heating  at  high  latitudes 
are  likely  to  be  more  variable  than  those  in  a  similar  experiment  at  the 
lower  latitudes.  We  discuss  the  results  in  the  context  of  the  ionospheric 
modification  (heating)  experiments  in  the  discussion  section.  The  theory 
of  parametric  excitation  of  E-region  modes  by  an  HF  pump  is  presented  in 
the  next  (second)  section,  the  lower  and  upper  E-region  modes  are 
discussed  in  the  third  and  fourth  sections.  The  results  are  discussed  in 
the  fifth  section  and  a  summary  is  given  in  the  sixth  section. 

2 .  THEORY 

We  consider  the  interaction  of  the  HF  pump  wave  with  the  low-frequency 
plasma  modes  (and  high  frequency  sidebands)  using  the  two-fluid  equations. 
The  coordinate  system  for  the  high-latitude  E-region  geometry  is  chosen 
so  that  the  earth's  magnetic  field  is  aligned  with  the  z-axis  (vertical), 
the  cross-field  Hall  current  (carried  by  an  equilibrium  electron  Hall 
drift,  v^)  is  along  the  y-axis  (east-west),  and  the  equilibrium  density 
gradient  is  along  the  x-axis  (north-south).  The  parallel  return  currents, 
JQ||,  are  represented  by  an  equilibrium  parallel  drift,  v^ | |  of  thermal 
electrons.  The  low  latitude  geometry  is  obtained  by  rotating  the  zx-plane 
by  90°  around  the  y-axis  so  that  the  z-axis  (horizontal)  points  towards 
geographical  north  and  the  x-axis  (vertical)  is  pointing  downwards.  The 
two-fluid  equations  are  [Braginskii,  1965] 
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In  the  above,  most  of  the  symbols  have  their  standard  meanings:  n,  m,  v, 

p,  T,  v  are  respectively  the  density,  mass,  velocity,  pressure  (scalar), 

temperature  (expressed  in  energy  units)  and  collision  frequency  of  the 

species  a  where  the  subscript  a  is  used  for  ions  (i)  and  electrons  (e). 

The  plasma  is  singly  ionized  (Z  =  1)  and  thus  the  equilibrium  electron 

plasma  density  (noe)  and  ion  plasma  density  ( nQ ^)  are  assumed  equal  and 

represented  by  nQ  (nQe  =  n^  =  nQ).  Q  is  the  heat  generation  term  and  q 

is  the  heat  flux  vector.  The  electron  collision  frequency  (vg)  in  the  E- 

region  is  dominantly  the  electron-neutral  collision  frequency  (v  )  as  the 

e-i  collisions  are  not  important  in  this  region  (vgn  >  Thus  in  the 

electron  energy  balance  equation,  the  dominant  contribution  to  the  heat 

generation  term  (Q  )  comes  from  the  e-n  collisions.  Further,  0  =  3m  /m 

°  e'  n  e  n 

v  is  the  inverse  of  the  energy  relaxation  time  for  the  elastic 
en 

collisions.  In  this  work,  we  do  not  consider  the  contributions  in  the 
electron  energy-balance  equation  coming  from  the  inelastic  collisions,  for 
simplicity,  though  these  could  be  of  importance  in  the  E-region  [Perkins 
and  Roble,  1978].  The  above  equations  are  supplemented  by  the  equations 
for  fields, 
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In  the  analysis  below,  we  have  neglected  the  temperature  effects  and 

the  ambient  magnetic  field  (Bq)  effects  on  the  high  frequency  modes.  Both 

of  these  approximations  are  made  for  simplicity  and  are  straightforward  to 

include  in  the  analysis.  As  a  consequence  of  these  approximations,  the 

high  frequency  modes  in  our  treatment  are  electron  plasma  oscillations, 

and,  the  pump  wave  frequency  (wq)  is  chosen  so  as  to  be  near  the  local 

electron  plasma  frequency  (wpe)*  1°  the  E-region  ionosphere,  the  electron 

plasma  frequency  is  larger  than  the  electron  cyclotron  frequency  (Q  ). 

5-3  e 

[For  normal  conditions,  for  n  >10  cm  ,  co  /S  >  2  while  during  the 

e  ~  pe  e  ~ 

periods  of  magnetic  activity,  this  ratio  may  be  larger  at  high  latitudes 
6—3 

e.g.,  the  ng  ~  10°  cm  ,  Wpg/^e  ~  6.]  Thus,  the  neglect  of  magnetic  field 

effects  on  HF  modes  is  not  a  good  approximation  under  the  normal  E-region 

conditions.  In  the  realistic  situation,  with  magnetic  field  effects 

included,  the  pump  wave  is  reflected  at  an  altitude  where  the  pump  wave 

2  2 

frequency  matches  the  local  upper-hybrid  frequency,  (where  =  w  g  + 

8^,  for  a  collisionless  cold  plasma),  rather  than  at  the  higher  altitude 

where,  w  ~  w  ,  as  has  been  assumed  in  our  treatment.  In  reality,  the 
o  pe 

reflection  occurs  even  at  a  somewhat  lower  altitude  due  to  dispersive 
effects.  A  systematic  study  of  the  inclusion  of  magnetic  field  effects  on 
the  parametric  instability  processes  has  been  carried  out  by  Weatherall  et 
al.  [1981].  Weatherall  et  al.  [1981]  have  shown  that  the  inclusion  of 
magnetic  field  effects  leads  to  nonlinear  frequency  shifts  that  affect  the 
unstable  wavenumber  spectra.  Our  analysis  below  focusses  primarily  on  the 
possibility  of  artificial  excitation  of  the  low  frequency  modes,  natural 
to  the  ambient  ionosphere,  in  the  presence  of  a  powerful  HF-pump  wave, 
and,  on  the  approximate  comparative  threshold  conditions  required  for  the 
excitations  of  these  modes  in  different  regions  of  the  E-region.  Thus, 
the  omission  of  the  above  effects  (Tq  =  0,  Bq  =  0  for  HF  modes)  may  not 
alter  importantly  the  conclusions  on  the  relative  threshold  estimates  and 
instability  conditions  for  the  various  modes  discussed  below.  We  intend 
to  include  these  effects  in  a  future  study.  In  the  reflection  region,  the 
electromagnetic  pump  wave  may  undergo  a  resonant  mode  conversion  into  an 
electrostatic  mode  [Piliya,  1967].  In  the  analysis  below,  we  have  assumed 
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an  electrostatic  pump.  The  pump  wave  (coq,  kQ)  is  given  by 
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where  nQ  is  the  undisturbed  ambient  ionospheric  plasma  density  and  n  is 
the  equilibrium  density  fluctuation  associated  with  the  HF  electrostatic 
pump  wave.  The  pump  wave  induced  oscillatory  electron  velocity  VeQ  is 
often  referred  to  as  the  "quiver  velocity".  The  interaction  of  the  pump 
with  the  low  frequency  perturbations  (co,  k)  generates  high-frequency 
sidebands  (co  ±  coq,  k  +  kQ,  with  co  <<  coq ) .  The  HF  sidebands  are  assumed 
electrostatic  and  are  described  by 
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where  subscript  H  refers  to  HF.  In  (6c),  the  terms  on  the  right  are  the 
nonlinear  terms.  The  low-frequency  modes  considered  in  this  paper  satisfy 
the  inequality,  co  <<  C0p^,  and  as  a  consequence  quasi-neutrality  is  assumed 
for  the  low-frequency  modes,  i.e., 
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The  ions  are  unaffected  by  the  HF  pump  and  are  described  by 
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The  magnetic  field  effects  on  ions  are  retained  in  the  above  as  they 
assume  importance  at  higher  altitudes.  The  ion  temperature  effects  have 
been  neglected  in  eqs.  (8)  for  simplicity  and  vid  is  the  equilibrium  ion 
drift  velocity.  The  electron  dynamics  on  the  slower  time  scale  is 
described  by 
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where  vt2=  (Tg/  me);  *L=  (eVTeo>;  D±=  (v^2/  Q2); 


2  2 


“aL=  (V  ~L*  ~ad)  ’  *e=  K7Te;  Kl=  (cX,Venevt  /SJe  );  (9f) 
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In  (9),  subscript  L  refers  to  the  low-frequency,  and,  it  has  been  assumed 

that  v  <<  2  and  v  >>  v  ..  Also  in  the  above,  q  is  the  heat  flow 

e  e  en  e i  - 

vector,  K  is  the  thermal  conductivity  and  c  is  the  coefficient  of  thermal 

A 

conductivity.  The  perpendicular  and  parallel  components  of  quantities 
with  respect  to  the  magnetic  field  are  denoted  by  subscripts  1  and  | | 
respectively.  The  nonlinear  heat  generation  term  [Q]NL  *n  (9d)  for 
electrons  arises  as  a  result  of  the  resistive  electron  heating  under  the 
action  of  high  frequency  fields  (m,n)  (eq.  9e)  where  (m,n)  represent  the 
interacting  modes.  Ue  are  considering  a  four-wave  process  in  which  the 
high-frequency  modes  are  the  pump  wave  and  the  two  side-bands  and  the 
fourth  mode  is  a  low  frequency  mode. 

The  standard  mode  coupling  analysis  is  followed  in  the  derivation  of 
the  dispersion  relation.  We  note  here  that  the  nonlinear  ponderomotive 
force  effects  have  been  neglected  in  (9)  as  they  are  usually  smaller  than 
the  nonlinear  thermal  coupling  effects  in  collisional  plasmas  if  the 
inequality,  \>e  >>  |w^|  ,  is  satisfied  [see  e.g.,  Papadopoulos  et  al., 

1983].  This  is  the  case  for  the  low-frequency  modes  considered  in  this 
work.  In  the  following,  we  present  results  of  the  stability  analysis  for 
different  modes. 

3.  LOVER  E-REGION  MODES 

We  first  consider  the  well  known  instabilities  associated  with  the 
equilibrium  Hall  current  (Lee  et  al.,  [1972]).  These  are  likely  to  occur 
predominantly  in  the  100-115  km  altitude  range.  In  this  region  magnetic 
field  effects  on  the  ion  dynamics  may  be  ignored  (v^n  »  SL). 
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(1)  Two-Stream  Instability  (TSI) 

This  instability,  also  known  as  the  Farley-Buneman  instability  (FBI), 
occurs  when  the  equilibrium  electron  Hall  drift,  v^,  exceeds  the  ion 
sound  velocity  [Farley,  1963;  Buneman,  1963].  This  instability  is 
discussed  in  the  auroral  and  equatorial  E-region  as  possibly  the  cause  of 
the  irregularities  associated  with  the  type  I  VHF  backscatter  radar  echoes 
[see  e.g.,  Schlegel  and  St. -Maurice,  1983;  Farley  and  Balsley,  1973].  Ue 
first  give  the  equation  describing  the  two-stream  instability  when  the 
ponderomotive  force  effects  were  included  but  the  nonlinear  Joule  heating 
effects  were  omitted  [Lee  et  al.,  1972],  i.e., 
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where 


(10a) 


From  (10),  it  is  readily  observed  that  the  pump  destabilizes  (stabilizes) 

the  two-stream  modes  when  6  <  0,  w  <w  (S  >  0,  &>>«).  The  pump 

o  pe  o  pe  _  v  r 

power  threshold  to  influence  the  two-stream  modes  is  ~  1  mW/ni  and  the 
pump  electric  field  threshold  is  -  0.2  V/m.  The  above  criterion  for 
instability  and  the  threshold  estimates  are  similar  to  the  case  of  the 
equatorial  electrojet  considered  by  Lee  et  al.  [1972].  The  high  latitude 
E-region  parameters  used  in  making  the  above  estimate  are,  v  -  4  x  1CT 
sec-1;  ~  2.5  x  10  see-  ;  v^  -  3.5  x  10  cm  sec-  ;  A  ~  3m.  As  has 

been  noted  before,  the  consideration  of  ponderomotive  force  nonlinearity 
and  neglect  of  the  thermal  nonlinearity  is  appropriate  for  situations  when 
Iw^l  >  vg.  In  the  case  of  the  two-stream  instability,  though,  the 
opposite  is  true  and  thus  the  dominant  nonlinear  coupling  is  provided  by 
the  nonlinear  electron  Joule  heating  effects,  in  which  case,  the  two- 
stream  instability  is  described  by 
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In  the  case  of  the  two-stream  modes,  0n  <  |  <*>e I  and  <0  ^  <  0.  The  pump 

destabilizes  (stabilizes)  the  two-stream  instability  for  5  <  0,  w  <  w 

o  pe 

(5  >  0,  o>o  >  w  ) .  The  threshold  pump  electric  field  (power)  to  influence 

the  two-stream  modes  when  the  nonlinear  Joule  heating  effects  are  dominant 

(for  the  above  mentioned  parameters)  is  -  60  mV/m  (-  5uU/m  ).  Ue  find 

that  the  threshold  values  of  pump  power  (electric  field)  to  influence  the 

two-stream  mode  are  significantly  lover  in  the  case  when  the  nonlinear 

Joule  heating  effects  are  dominant  (eq.  11)  than  the  case  when  the 

nonlinear  ponderomotive  force  effect  are  important  (eq.  10).  The  scale 

sizes  generated  by  the  two-stream  instability  approach  -  a  few  meters 

(which  is  close  to  the  ion  mean  free  path,  <  meter).  For  these 

wavelengths,  the  ion  kinetic  equation  should  be  used  to  include  ion 

Landau-damping  effects.  The  fluid  equation  approach  appears  to  give 

reliable  estimates  for  wavelengths  larger  than  (for  X^>  2m).  The  above 

discussion  also  holds  for  most  part  to  the  low-latitude  situations.  At 

low  latitudes,  the  electrojet  peak  altitude  is  -  105km  (in  contrast  to  the 

high  latitude  situation  where  it  is  -  110km).  The  peak  value  of  electron 

drift  at  low  latitudes  generally  exceeds  the  local  ion-sound  speed 

nominally  but  not  too  greatly  (c  <  <  2c  );  whereas,  at  high 

peak  s  ai  s 

latitudes,  v^  may  at  times  exceed  cg  several  limes  over. 

We  have  shown  above  that  in  the  presence  of  an  HF  pump  wave  of 
2 

strength  5yW/m  (pump  electric  field  60  mV/m),  the  two-stream  modes  in  the 
E-region  (at  altitudes  near  110  km  for  high-latitudes  and  near  105  km  for 
low  latitudes)  may  be  excited  for  sub-threshold  Hall  drifts  (v^  <  cg)  in 
the  overdense  region  («Q  <  wpe^‘  ^  was  f°unc*  t^iat  t^ie  nonlinear  Joule 
heating  is  the  dominant  coupling  mechanism  for  these  modes  in  the  E-region 
ionosphere  (as  compared  with  the  ponderomotive  force  nonlinearity). 
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(2)  Gradient  Drift  Instability  (GDI) 

This  instability  occurs  when  the  electron  Hall  drifts  are  sub-sonic 
(v^  <  cg)  in  the  presence  of  an  equilibrium  transverse  density  gradient 
(which  for  instability  should  be  in  the  same  sense  as  the  primary 
polarization  electric  field  responsible  for  the  electron  Hall  drift) 
[Simon,  1963;  Hoh,  1963].  The  gradient-drift  instability  favors 
relatively  longer  scale  size  irregularities,  ~  100's  meters  in  the  E- 
region  [Rogister  and  D'Angeleo,  1970].  A  two-step  theory  proposed  by 
Sudan  et  al.  [1973]  for  the  equatorial  case  has  been  frequently  invoked  to 
interpret  the  type  II  E-region  VHF  backseat  ter  echoes  at  the  high 
latitudes  also  [Greenwald,  1974;  Schlegel  and  St. -Maurice,  1983].  In  the 
two-step  theory,  the  large  scale  size  irregularities  generated  by  the 
gradient-drift  instability  grow  to  large  amplitudes  and  then  nonlinearly 
generate  secondary  smaller  scale  size  (~  3  m's)  via  the  two-stream  and/or 
gradient-drift  instability  mechanisms.  Sharma  [1989]  has  studied  the 
gradient-drift  instability  in  presence  of  an  HF  pump  in  the  equatorial  E- 
region  case  including  the  thermal  nonlinearity.  We  present  results  for 
the  high-latitude  case  below.  However,  first  we  present  results  of  the 
calculation  for  the  situation  when  the  dominant  nonlinearity  is  due  to  the 
ponderomotive  force  effects  (weakly  collisional  case,  vg  <  lo^l)  and  when 
the  nonlinear  ohmic  heating  effects  are  neglected  (the  case  considered  by 
Lee  et  al.  [1972]).  The  dispersion  relation  for  the  gradient  drift 
instability  in  this  case  is  given  by 


«(1  +  a)  -  k1*vd  +  iDj_k 2  *  i  v.nko« 


2v  | v  |2k2w2  6 
e1  o'  Ip 
_ *e 

n2,,2  2  2. 

2(6  +  v  (o  ) 

e  e  o 


(12') 


where  k  =  e  • ( e  x  k.)/k, 
o  z  ~n  l 


From  (12'),  one  may  readily  calculate  the  real  frequency,  co^,  and  growth 
rate,  |y|»  by  substituting  «  =  «r  +  i | t I  ( I T I  <  «r).  This  expressions  is 
similar  to  the  one  obtained  by  Lee  et  al.  [1972]  and  yields  similar 
criteria  for  influencing  the  gradient-drift  instability  by  the  external 
pump  wave,  i.e.,  destabilization  (stabilization)  occurs  when  6  <  0,  a>o  < 
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4  -1 

«  (5  >  0,  to  >  (u  ).  For  the  low  latitude  parameters,  v  ~  4  x  10  s 

pe  o~  _jPe  4  _1  e 

~  2.5  x  10  s  ,  ~  10  cm  s  ,  and  ~  100m,  the  pump  electric  field 

required  for  influencing  the  gradient  drift  mode  in  the  E-region  is  -  0.2 

V/m,  a  value  on  the  same  order  of  magnitude  as  that  obtained  by  Lee  et  al. 

[1972]. 

For  the  collisional  case  (vg  >  |aijJ)>  appropriate  for  the  E-region 
ionosphere,  we  need  to  include  the  ohmic  heating  nonlinearity.  The 
dispersion  relation  in  this  case  may  be  obtained  from 


....  .  2  .  vt  kl  IV1(Vn)1<ko-lve> 

“L(l+a)  -  *  <Dlkl  *  yL>  *  - ~  —  - 


(12) 


From  (12)  the  growth  rate  for  the  gradient-drift  instability  in  the 
presence  of  an  HF  pump  wave  may  be  written  as 
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From  (12),  we  find  that  the  pump  is  destabilizing  (stabilizing)  when  S  < 
0,  w  <«  (6  >  0,  w  >w  )  (note  that  for  these  modes  «  .  <  0  and  for 

the  linearly  unstable  situation,  k  <  0).  The  threshold  pump  field 

4  -1  0  4-1 

requirement  for  v, .  -  1  x  10  cm  s  ,  X,  -  100m,  v  -  4  x  10  s  ,  v. 

«  I  qJL  J.  6  in 

2.5  x  10J  s'  ,  etc.,  is  15  mV/m.  This  estimate,  for  the  parameters 
mentioned,  is  applicable  to  both  the  high  and  low  latitudes.  Ve  find  that 
the  threshold  pump  electric  field  in  the  E-region  to  influence  the 

gradient  drift  mode  is  easily  down  by  an  order  of  magnitude  (or  more)  in 
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the  case  when  the  nonlinear  ohmic  heating  effects  were  included  in  the 

analysis  (eq.  12),  compared  with  the  case  when  collisional  effects  were 

neglected  and  the  nonlinearity  was  due  to  the  ponderomotive  force  effects 

(eq.  12').  The  stability  theory  here  has  used  the  local  approximation  in 

deriving  the  dispersion  relation.  This  approach  is  valid  when  »  1. 

For  scale  sizes  such  that  k^L^  <  1,  a  nonlocal  approach  should  be 

considered  for  the  analysis.  In  the  appendix,  ve  consider  the  effects 

arising  from  such  a  treatment.  For  the  medium  scale  size  modes  such  that 

2 

k^L^j  >10  (X  ~  100  m)  the  local  theory  estimates  may  be  considered 

representative. 

4.  UPPER  E-REGION  MODES 

At  higher  E-region  altitudes  at  high  latitudes,  the  effect  of  parallel 
currents  introduces  the  possibility  of  additional  plasma  modes.  These 
modes  are  likely  to  be  favored  at  altitudes  >  110  km.  In  this  region 

decreasing  ion-neutral  collision  frequency  implies  inclusion  of  ion 
magnetization  effects  (y^  <  52^).  But,  depending  on  the  ambient 
conditions,  the  regions  of  occurrence  for  the  lower  E-region  modes  (two- 
stream  instability  and  gradient-drift  instability)  and  the  upper  E-region 
modes  may  often  overlap. 

(1)  Ion  Acoustic  (IA)  Mode 

At  altitudes  above  >  110  km,  nearly  transverse  propagating  ion- 

acoustic  modes  may  be  driven  unstable  by  a  combination  of  parallel  and 
perpendicular  equilibrium  electron  drifts  [Kaw,  1973;  Chaturvedi  et  al., 
1987b;  Villian  et  al.,  1987].  There  are  indications  that  some  auroral 
backscatter  radar  observations  may  be  interpreted  as  the  ion-acoustic 
modes  [Villian  et  al.,  1987;  Haldoupis  et  al.,  1986].  The  dispersion 
relation  in  this  case  may  be  obtained  from 


k2c  2  k2c  2h 
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where  [D:k2]  s  (D  k  2+D, ,k. ,2); 
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The  criterion  for  instability  in  this  case  is  as  follows.  For  modes  with 
<  0  (the  minimum  condition  for  linear  ion-acoustic  instability),  the 
pump  is  further  destabilizing  (stabilizing)  for  5  <  0,  (0Q  <  w  (5  >  0, 

>  « ).  For  modes  with  «  .  >0  (linearly  stable  case),  the  reverse  is 

pe  eL  J  2  1 

true.  The  threshold  pump  field  in  this  case,  for  \>.  -  3.6  x  l(r  s  ,  v  - 

3-1  7  -1  e 

5.3  x  10  s  ,  Tg  ~  O.leV,  ~  1.3  x  10  cm  s  ,  etc.,  is  30  mV/m.  A 

more  precise  determination  of  the  threshold  for  the  HF  induced  ion- 

acoustic  instability  may  be  obtained  by  using  the  kinetic  ion  equations  so 

as  to  include  the  ion-Landau  damping  effects.  The  above  estimate  is  a 

good  approximation  for  scale  sizes  larger  than  the  ion  mean  free  path  (^  1 

meter) . 

(2)  Electrostatic  Ion  Cyclotron  (EIC)  Mode 

There  is  considerable  speculation  that  at  the  upper  E-region  altitudes 
(>  130  km),  the  generation  of  EIC  modes  by  parallel  currents  may  result 

in  yet  another  type  of  auroral  backscatter  radar  return  [D'Angelo,  1973; 
Chaturvedi,  1976;  Fejer  et  al.,  1984;  Satyanarayana  et  al.,  1985; 
Providakes  et  al.,  1985;  Prikryl  et  al.,  1988;  Villain  et  al.,  1987]. 
This  instability  is  favored  at  altitudes  where  vi  <<  Cj.  Chaturvedi  et 
al.  [1987a]  and  Laxmi  and  Tripathi  [1988]  have  considered  the  possibility 
of  HF-induced  ion-cyclotron  mode  excitation  but  both  these  papers  used  the 
ponderomotive  force  nonlinearity.  The  nonlinear  Joule  heating  of 
electrons  may  be  the  dominant  coupling  mechanism  in  this  case  also,  as 

>  | (o^|  is  likely  to  be  satisfied.  The  dispersion  relation  may  be  obtained 
from 
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It  may  be  noted  that  the  form  of  pump-associated  terms  in  the  case  of  the 

EIC  mode  (eq.  14)  is  similar  to  the  one  in  the  ion-acoustic  mode  case  (eq. 

13).  This  is  not  surprising  as  the  electron  dynamics  in  the  two  cases  are 

quite  similar.  The  instability  criterion  for  the  HF-induced  EIC  mode  is 

thus  similar  to  the  ion-acoustic  mode  case.  But  the  occurrence  of  the  EIC 

mode  is  likely  at  higher  altitudes  when  ^  »  \^;  whereas,  the  ion 

acoustic  mode  may  be  excited  at  lower  altitudes,  such  that  8.  <  v.  .  The 

l  -  in 

threshold  pump  field  for  the  EIC  instability,  for  v^n  -  10  s'  ,  is  - 
5mV/m.  Due  to  the  neglect  of  ion  temperature  and  the  use  of  fluid 
equations  for  simplicity,  the  theory  presented  here  has  not  considered  the 
case  of  higher  harmonics  («n  -  nQ^).  The  use  of  ion  kinetic  equations 

ensures  inclusion  of  these  effects  and  ion-Landau  damping  effects 
[Chaturvedi  et  al.,  1987a],  and  is  planned  for  future  work. 


(3)  Current  Convective  Instability  (CCI) 

This  instability,  which  occurs  in  the  presence  of  a  parallel  current 
and  a  transverse  density  gradient,  has  been  applied  to  the  generation  of 
high-latitude  F-region  irregularities  in  the  intermediate  scale  size  range 
(~  100  m's)  [Ossakow  and  Chaturvedi,  1979).  In  the  lower  E-region  (100  - 
115  km),  the  growth  rate  of  the  gradient  drift  instability  is  much  larger 
than  the  current-convective  instability  [Chaturvedi  and  Ossakow,  1981], 
but  at  higher  altitudes  in  the  E-region  120  km)  the  current-convective 
instability  may  occur  and  generate  longer  scale  size  (  100  m's) 
irregularities.  The  dispersion  relation  of  the  current-convective 
instability  may  be  obtained  from 
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In  absence  of  the  pump  wave,  eq.  (15)  yields  the  familiar  dispersion 

relation  for  the  current-convective  instability  [Ossakow  and  Chaturvedi, 

1979].  In  this  case,  ©n  >  jw^|  and  the  criterion  for  HF  induced 

destabilization  (stabilization)  for  the  current-convective  instability  is 

S  <  0,  to  <  to  (5  >  0,  0)  >w  ).  The  threshold  pump  field  required  in 

o  pe  o  pe 
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this  case,  for  v  /Q  ~  10 

5  e  e  5 

l(r  cm  s  ,  Ljj  ~  5  x  10’’ 


"4  /k.  -  (m  /m.)1/2,  v./Q.  -  10~3, 

i  v  e  '  l  i  1 


m,  is  ~  2mV/m. 
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This  instability  is  similar  to 
the  gradient  drift  instability  and  generates  relatively  longer  scale  size 
irregularities.  The  above  analysis  has  been  carried  out  using  the  local 
approximation  which  is  valid  for  wavelengths  satisfying,  k^L^  »  1.  At 
longer  wavelengths,  such  that  k^L^  ^  1,  non-local  effects  are  important 
and  have  been  considered  in  the  appendix  for  the  case  of  the  gradient 
drift  instability.  For  the  scale  sizes  of  interest  (A^  -  100m,  for  ~  5 
km,  k^Ljj  -  102),  the  local  theory  expression  may  be  considered  adequate. 


5.  DISCUSSION 

In  the  above,  we  have  described  several  low  frequency  plasma 
instabilities  which  may  be  influenced  by  the  presence  of  a  large  amplitude 
HF  pump  in  the  high  and  low  latitude  E-region.  There  have  been  several 
artificial  ionospheric  modification  (AIM)  heating  experiments  resulting  in 
the  generation  of  small  scale  size  plasma  irregularities  in  the  high 
latitude  E-region.  The  various  experimental  observations  are  indicative 
of  apparently  a  variety  of  processes  occurring  in  this  region  during  the 
ionospheric  heating  experiments  [Stubbe  and  Kopka,  1983;  Hibberd  et  al., 
1983;  Coster  et  al.,  1985;  Djuth  et  al.,  1985;  Hoeg  et  al.,  1986],  Some 
features  of  these  observations  have  been  explained  in  terms  of  the  well 
known  plasma  instability  mechanisms  associated  with  the  HF  pump  wave  such 
as  the  resonance  instability  [Vaskov  and  Gurevich,  1977;  Inhester,  1982], 
and  the  thermal  parametric  instability  [Grach  et  al.,  1977;  Das  and  Fejer, 
1979;  Dysthe  et  al.,  1983].  Some  aspects  of  the  observations  point 
towards  a  role  that  the  natural  plasma  modes  of  this  region  may  be 
playing.  For  instance,  the  observed  radar  backscatter  spectra  of  the 
artificially-generated  irregularities  (Al)  show  similarities  to  that 
obtained  in  the  case  of  natural  irregularities,  and,  in  some  observations 
the  effects  seen  associated  with  the  artificially  excited  irregularities 
are  similar  to  those  in  the  case  of  the  natural  irregularities  (e.g., 
enhanced  heating).  In  the  following,  we  discuss  the  E-region 
instabilities  discussed  in  the  previous  sections  in  view  of  the  previous 
and  possible  future  ionospheric  modification  (heating)  experiments. 
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(a)  Two  Stream  Instability  -  This  instability  has  been  discussed  in  the 
context  of  ionospheric  heating  experiments  at  high  latitudes  by  Ribberd  et 
al.  [1983]  and  Hoeg  et  al.  [1986].  The  linear  instability,  without  the 
pump  wave,  requires  that  the  ambient  d.c.  electric  field  be  larger  than  - 
20  mV/m  (i.e.,  the  relative  electron  Hall  drift  be  larger  than  the  ion 
sound  speed,  -  350  m/s).  The  Hall  current  peak  is  around  ~  110  km 

altitude  in  the  high-latitude  E-region  and  this  instability  is  favored 
around  the  region  of  peak  current  intensity.  At  higher  altitudes,  the 
relative  drift  between  electrons  and  ions  decreases  as  the  ion 
magnetization  effects  (and  the  ion  drift  effects)  increasingly  assume 
importance.  In  regions  where  the  ionosphere  is  stable  to  the  linear  two 
stream  instability, ,  a  contribution  from  another  source  of  free  energy  may 
create  conditions  for  the  onset  of  the  instability.  In  an  ionospheric 
heating  experiment,  an  HP  pump  wave  may  thus  excite  the  two  stream 
instability  in  the  regions  even  where  the  current  is  linearly  sub¬ 
threshold.  In  Hoeg  et  al.  [1986],  one  set  of  modification  experiment 
observations  at  higher  E-region  altitudes  is  described  as  two-stream-like 
at  a  time  when  the  natural  two  stream  instability  was  observed  at  lower 
altitudes.  The  threshold  pump  electric  field  required  for  the  two  stream 
instability  is  -  60  mV/m  which  is  on  the  order  of  the  pump  wave  strength 
in  these  experiments  in  the  E-region  (^90  mV/m).  The  natural  type  I 
irregularities  have  often  been  associated  with  enhanced  electron 
temperatures  in  the  high  latitude  E- region  [Schlegel  and  St. -Maurice, 

1981].  In  the  presence  of  the  HF  pump  wave,  the  electron  temperature  is 

2 

enhanced  by  an  amount  proportional  to  (m  V  /3)  where  m  is  the  mass  of 

n  eo  n 

neutrals  and  VeQ  is  the  electron  quiver  velocity.  An  artificial 
ionospheric  modification  (heating)  experiment  may  be  carried  out  to 
reproduce  the  phenomenon  of  increased  Tg  associated  with  the  type  I 
natural  high-latitude  E-region  irregularities  over  the  above  mentioned 
enhancement.  The  growth  times  of  the  two-stream  instability  generated 
irregularities  are  on  the  order  of  -0.1  sec. 

Furthermore  there  appear  to  have  been  no  heating  experiments  at  the 
low  latitudes.  An  artificial  excitation  of  type  I  modes  during  sub¬ 
threshold  conditions  and/or  at  altitudes  away  from  the  peak  electrojet 
layer  may  aid  in  testing  the  two-stream  instability  interpretation  of  the 
natural  type  I  echoes. 
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(b)  Gradient  Drift  and  Current-Convective  Instabilities  -  The  gradient 

drift  instability  may  occur  when  the  d.c.  electric  field  is  sub-threshold 

for  the  two-stream  instability  (<  20  mV/m),  requires  a  transverse  density 

inhomogeneity,  and  has  a  growth  time  period  on  the  order  of  ~  tens  of 

secs.  Hibberd  et  al.  [1983]  discarded  the  two-stream  instability 

interpretation  of  their  observations  on  account  of  the  slow  time  period 

associated  with  the  irregularities  in  their  experiment  ( ~  30  sec),  which 

was  much  larger  than  the  two-stream  instability — associated  time  scale  of 

-  0.1  sec.  Later,  for  some  other  ionospheric  modification  experiments, 

the  gradient  drift  instability  had  been  considered  as  a  possible  mechanism 

[Djuth  et  al.,  1985;  Hoeg  et  al.,  1986].  This  instability  linearly 

generates  larger  scale  sizes  (~  100's  m)  which  may  grow  to  larger 

amplitudes  and  then  nonlinearly  generate  the  smaller  scales  via  the  two- 

step  process  [Sudan  et  al.,  1973].  The  pump  wave  threshold  electric  field 

required  for  the  gradient-drift  instability  in  the  high-latitude  E-region 

is  ~  15  mV/m  for  an  altitude  ~  110  km  and  a  d.c.  electric  field  strength 

of  -  10  mV/m.  Thus  a  combination  of  ambient  d.c.  electric  field  and  pump 

wave  of  less  than  the  peak  strengths  reported  in  the  above  mentioned 

modification  heating  experiments  seems  capable  of  readily  exciting  the 

gradient  drift  instability.  At  higher  altitudes,  the  electrojet  strength 

ebbs  and  larger  pump  wave  strength  may  be  required  for  the  gradient  drift 

mode  excitation.  But  at  high-latitudes,  the  presence  of  field-aligned 

currents  (FAC)  may  result  in  a  pump  wave-induced  gradient  instability,  the 

current-convective  instability  (CCI),  even  at  higher  altitudes  (outside 

the  electrojet  region).  The  pump  wave  threshold  for  the  current 

convective  instability  at  ~  150  km  is  >  2  mV/m  where  a  presence  of  field 

2 

aligned  current  of  strength  -  10  yA/m  has  been  assumed.  The  current 

convective  instability  generates  longer  scale  size  irregularities  (~  100's 
m)  and  the  associated  growth  periods  are  ~  10 's  secs.  A  cascade  to 

smaller  scale  sizes  via  the  two-step  mechanism  may  then  result  in  meter- 
scale  irregularities  detectable  by  a  VHF  radar.  Thus  in  an  ionospheric 
modification  (heating)  experiment,  it  may  be  possible  to  observe  HF  pump 
wave  induced  irregularities  well  above  the  electrojet  region  and  still 
having  slow  growth  time  scales.  The  spectra  of  irregularities  thus 
generated  will  be  of  type  II  if  detected  by  a  VHF  backscatter  radar. 
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A  similar  experiment  may  be  carried  out  at  low  latitudes  to 
artificially  excite  (or  suppress)  the  type  II  irregularities  at  low 
latitudes.  A  generation  of  type  II-like  irregularities  in  the  regions 
away  from  the  electrojet  layer  may  thus  aid  in  testing  the  model  of 
natural  type  II  echoes  in  terms  of  the  gradient  drift  instability  and  the 
two-step  theory. 

For  the  gradient  drift  instabilities,  an  important  consideration  at 
large  wavelengths  is  the  modification  of  the  growth  rate  and  the  threshold 
criterion  due  to  the  non-local  effects.  The  local  theory  approximation  is 
generally  valid  when  k^L^  >>  1.  In  the  E-region,  the  gradient  scale 
length,  I^j  ~  6  kms  and  for  scale  sizes  of  interest  (X^  ~  100's  m),  k^L^  ~ 
10^,  and  thus  the  nonlocal  effects  may  still  be  regarded  as  small.  For 
smaller  values  of  L^,  and/or  for  larger  X's  such  that  k^LN  ~  1,  a  nonlocal 
analysis  is  required.  In  this  case  (k^L^  ~  1),  the  growth  rates  are 
reduced  and  larger  pump  thresholds  are  required. 

(c)  Ion  Acoustic  and  Electrostatic  Ion  Cyclotron  Modes  -  The  observations 
of  natural  plasma  irregularities  in  the  high  latitude  E-region  ionosphere 
have  indicated  a  presence  of  apparently  many  types  of  plasma 
irregularities  in  contrast  to  the  equatorial  case  where  most  of  the 
observations  have  been  interpreted  in  terms  of  the  two-stream  modes  (type 
I  radar  spectra)  and  gradient  drift  modes  (type  II  spectra).  The 
irregularity  classification  of  high  latitude  observations  is  still 
emerging,  based  on  the  VHF  and  HF  backseat  ter  radar  spectra  and  in-si tu 
rocket  observations.  The  presence  of  low-frequency  ion  modes  in 
association  with  the  parallel  currents  at  high  latitudes  has  been 
considered  as  a  possible  interpretation  of  many  E-region  radar 
observations.  For  example,  the  EIC  mode  has  been  invoked  by  Fejer  et  al. 
[1984],  Prikryl  et  al.  [1988]  and  Villain  et  al.  [1987]  for  some  of  their 
observations.  Villain  et  al.  [1987]  have  mentioned  the  ion  acoustic  mode 
for  some  of  their  observations.  Thus  it  is  expected  that  the  ionospheric 
modification  experiments  at  high  latitudes  will  provide  a  richer  variety 
of  observed  phenomena  than  expected  at  the  lower  altitudes. 

In  regions  of  the  upper  E-region  electrojet  layer,  as  mentioned 
earlier,  the  threshold  condition  for  the  excitation  of  the  two  stream 
instability  may  not  be  met  because  of  ion  magnetization  effects.  It  has 
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been  suggested  that  an  obliquely  propagating  ion  acoustic  wave  (with  a 
finite  parallel  wave  vector)  will  be  unstable  to  a  combination  of 
transverse  and  parallel  electron  drifts  [Chaturvedi  et  al.,  1987b;  Villain 
et  al.,  1987]  and  at  higher  altitudes  the  field-aligned  currents  will 
destabilize  EIC  modes  [D'Angelo,  1973;  Chaturvedi,  1976;  Satyanarayana  et 
al.,  1985;  Providakes  et  al.,  1985].  In  the  ionospheric  modification 
(heating)  experiments,  an  artificial  excitation  of  these  modes  is  possible 
by  means  of  the  HF  pump  wave.  For  the  upper  electrojet  regions  where  ~ 
Q..,  the  EIC  mode  is  heavily  damped  due  to  ion-neutral  collisional  effects, 
but  the  excitation  of  the  ion  acoustic  mode  is  possible.  The  threshold 
pump  electric  field  required  for  the  ion  acoustic  instability  at  -  120  km 
altitude  is  ~  30  mV/m  which  is  in  the  range  of  the  artificial  ionospheric 
modification  (AIM)  heating  experiments  mentioned  above.  At  higher 
altitudes,  the  EIC  mode  could  be  excited  artificially  for  similar 
threshold  values.  The  EIC  generation  is  also  possible  at  F-region 
altitudes  [Chaturvedi  et  al.,  1987a]. 

Both  the  ion  acoustic  and  EIC  modes  are  obliquely  propagating  modes. 
The  natural  E-region  irregularity  radar  observations  at  high  latitudes 
have  often  detected  returns  at  large  aspect  angles  implying  a  finite 
parallel  wavenumber  associated  with  the  irregularities.  A  signature  of 
these  modes  would  be  a  similar  aspect  angle  dependence  observed  by  radar 
detected  artificial  generated  irregularities  during  ionospheric  heating 
experiments  in  the  upper  E-  and  F-regions.  Further,  at  higher  altitudes, 
the  nonlinearly  developed  EIC  turbulence  is  expected  to  result  in 
energetic  ions  (at  150  -  200  km)  or  possibly  ions  with  non-thermal 

distributions  (above  250  kms)  [Satyanarayana  et  al.,  1989].  Detection  of 
super-thermal  ions  could  be  used  as  a  possible  diagnostic  for  the 
artificially  induced  electrostatic  ion-cyclotron  wave  turbulence. 

The  generation  of  ion  acoustic  and  EIC  modes  in  the  presence  of  a 

field  aligned  current  is  possible  in  both  the  overdense  and  underdense 

regions,  depending  on  the  altitude  of  the  resonance  layer.  In  the 

linearly  stable  situation,  the  underdense  region  is  favored  for  the  pump 

wave-induced  instability  (wq  >  w  ).  But  in  presence  of  a  field  aligned 

current,  for  the  marginally  stable  or  unstable  conditions,  an  HF  pump  wave 

is  further  destabilizing  to  the  ion  acoustic  and  EIC  instabilities  in  the 

overdense  region  («  <  u  ).  The  time  scale  of  growth  for  these  modes  is 

o  pe 
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<  1  sec.  The  mode  excitation  by  this  mechanism  in  the  underdense  and 
overdense  regions  may  have  some  relevance  to  the  observations  of  Hoeg  et 
al.  11986]. 

6.  SUMMARY 

In  conclusion,  we  have  found  that  by  using  a  high  frequency  pump  wave 
at  high  latitude  E  region  altitudes,  the  following  observable  effects  may 
be  induced  in  an  artificial  ionospheric  modification  experiment. 

(i)  For  a  situation  when  the  electrojet  strength  is  sub-threshold  (v^ 

<  c  ),  or,  in  a  situation  when  at  the  peak  electrojet  altitudes  (at  -  110 

s 

km)  the  two  stream  instability  may  be  operative  (v^  >  c  )  but  at  the  same 
time  at  higher  altitudes  the  relative  electron  drift  may  be  sub-threshold 
due  to  the  ion  magnetization  effects  (lve^  -  I  <  cg),  a  pump  wave 
electric  field  strength  of  -  60  mV/m  may  excite  small  meter  scale  size 
irregularities,  observable  by  a  VHF  backseat  ter  radar  (type  I  spectra)  on 
a  time  scale  of  ~  0.1  sec.  These  irregularities  are  expected  to  result  in 
enhanced  electron  temperatures  in  the  region. 

(ii)  The  pump  wave  electric  field  of  strengths  ~  10  mV/m  could 
generate  medium  scale  (~  100's  meters)  irregularities  on  the  time  scales 
of  ~  tens  of  seconds  at  the  lower  and  upper  E-region  altitudes  via  the 
gradient  drift  and/or  current-convective  instabilities.  These 
irregularities  may  then  nonlinearly  generate  smaller  scale  sizes  (~ 
meters)  via  the  two  step  mechanism.  The  small  scale  size  irregularities 
are  detectable  via  a  VHF  backscatter  radar  (type  II  spectra).  The  primary 
irregularities  are  also  expected  to  result  in  the  scintillations  of  RF 
signals. 

(iii)  The  pump  wave  electric  field  strengths  of  ~  10  -  30  mV/m  could 
directly  generate  small  scale  size  (~  meters)  irregularities  via  the  ion- 
acoustic  instability  in  the  altitude  region  110-130  kms  and  via  the 
electrostatic  ion-cyclotron  instability  at  altitudes  *  130  km.  The 
detection  of  these  irregularities  is  possible  by  a  VHF  or  HF  backscatter 
radar.  These  modes  have  a  finite  parallel  wave  vector  and  a  study  of 
aspect  angle  dependence  of  radar  returns  could  be  used  to  reveal  this 
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feature.  At  higher  altitudes,  as  a  consequence  of  the  nonlinear 
development  of  the  EIC  instability,  the  detection  of  energetic  ions  would 
be  a  signature  of  these  modes. 

(iv)  Similar  experiments  may  be  carried  out  at  low  latitudes  also. 
These  could  include  generation  of  irregularities  (types  I  and  II)  in 
regions  away  from  where  the  peak  electrojet  strength  occurs  (near  105km 
altitude).  At  altitudes  separated  from  the  peak  electrojet  layer,  or, 
during  sub-threshold  electrojet  strengths,  an  artificial  excitation  of 
type  I  irregularities  will  be  a  test  of  the  two-stream  instability  model. 
Similarly,  a  generation  of  type  II  irregularities  at  altitudes  separated 
from  the  main  electrojet  region  will  verify  the  gradient-drift  instability 
and  two-step  theory  models  for  equatorial  electrojet  irregularities.  The 
strength  of  the  pump  electric  field  required  in  these  cases  is 
approximately  60mV/m  and  15mV/m  for  generation  of  the  type  I-like  and  type 
II-like  irregularities,  respectively. 
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APPENDIX 


NONLOCAL  E  x  B  INSTABILITY  IN  THE  E-REGION 

In  the  E-region,  there  have  been  observations  of  long-wavelength 
irregularities,  which  have  been  associated  with  the  electrojet  Hall 
currents  at  both  the  equator  [Kudeki  et  al.  1982,  Pfaff  et  al.  1987]  and 
at  high-latitudes  [Basu,  1988].  In  this  region,  the  observed  wavelengths 
are  ^  1  km  while  the  driving  density  gradient  scale  size  is  in  the  range 
of  -  6-10  kms.  Thus,  typically,  in  the  E-region,  for  the  long 
wavelengths,  kL^  -  10.  It  is  therefore  of  interest  to  examine  this 
instability  with  the  nonlocal  effects  included  at  long-wavelengths  in  the 
presence  of  the  pump  wave. 

In  this  section,  we  present  a  simple  model  calculation  for  the 
parametric  instability  of  the  GDI  by  an  HF-pump  to  illustrate  the 
importance  of  nonlocal  effects  on  this  instability.  We  consider  a 
simplified  density  model 


n  (x)  =  J 
ov  n. 


x  <  0 
x  >  0 


(Al) 


As  before  the  pump  and  the  high  frequency  modes  are  assumed  to  be 
electrostatic.  The  simplified  ion  and  electron  equations  are  given  as 
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We  have  neglected  the  diffusion  (and  recombination)  damping  effects  for 
simplicity  here.  We  use  eq.  (9)  in  (A3)  to  arrive  at 
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We  have  used  the  local  approximation  for  the  high-frequency  side-bands  in 
obtaining  (A4).  This  has  been  done  for  analytical  simplicity  and  may  be  a 
good  approximation  in  the  case  of  an  ES  pump  such  that  kQ  >  k^.  This 
enables  one  to  obtain  a  differential  equation  from  (A2)  and  (A4) 
which  describes  the  nonlocal  E  x  B  gradient-drift  mode, 
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We  solve  eq.  (A6)  by  using  the  method  described  by  Chandrasekhar  [1961] 
for  the  Rayleigh-Taylor  instability  of  a  heavy  fluid  resting  on  top  of  a 
lighter  fluid  against  acceleration.  A  similar  analysis  for  the  E  x  B 
instability  in  F-region  has  been  carried  out  by  Huba  and  Zalesak  f 1983 ] 
without  the  HF  pump  wave  effects.  A  solution  for  the  potential  $  is 
assumed  of  form 


♦  =  ^e  y  +  $2e  y  (A7) 

For  solutions  satisfying  the  boundary  conditions,  <f>  ->  0  as  |x|  ->  +  ®,  we 

get 
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x  >  0 

Moreover,  at  x  =  0,  the  potential  is  continuous,  thus 


(A8) 


=  *2  (A9) 

Integrating  Eq.  (A6)  across  the  boundary  (x  =  0)  and  making  use  of  (A7)  - 
(A9),  ve  obtain  the  nonlocal  dispersion  relation  for  the  GDI  in  the  E- 
region  in  the  presence  of  a  pump  wave 
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For  the  case,  ©n  >  I WjJ  >  ve  find  that  the  nonlocal  growth  rate  is 
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where  subscript  m  (=  1,2)  refers  to  the  magnitude  of  quantities  in  the  two 
regions. 

We  note  from  (All)  that  the  pump  is  destabilizing  if  ^  <  C  and  (J^n^l 
>  | ^2n2 | ,  which  would  roughly  correspond  to  the  earlier  obtained  criteria 
of  S  <  0  in  the  local  approximation.  We  note  that  ($m)  is  a  small  number 
in  only  one  region  (because  is  chosen  to  be  resonant  only  in  one 
region,  say,  region  1,  with  «Q  =  wpl^*  Then,  the  pump  is  no  longer 
resonant  in  the  region  2  and  therefore  |  I  >>  |  |  .  Thus,  in  this 
situation,  >  I ^n2 ^  *s  sat^s^ie^’  The  local  theory  expression  for 
the  growth  terms  for  the  GDI  in  the  presence  of  an  HF  pump  wave  (eq.  12a) 
may  be  obtained  from  (All)  by  substituting  for  ky((n^  -  n2)/(n^  +  n2)) 
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and  using  the  approximation  n^  >  n£.  Another  consequence  of  the  nonlocal 
theory  is  that  at  long  wavelengths,  the  growth  rate  is  proportional  to  the 
wavenumber  and  thus  the  pump-induced  contributions  are  smaller  at 
larger  scale  sizes.  Therefore,  larger  pump  wave  powers  are  required  to 
achieve  effects  comparable  to  the  case  of  small  scale  sizes.  Thus,  it  is 
more  difficult  to  influence  the  large-scale  size  interchange  process  than 
it  is  to  influence  the  small-scale  sizes  by  an  external  pump  wave. 
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